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Neutrino-nucleus reactions in the energy range 1-100 MeV 
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We review some salient aspects of calculations of the neutrino-nucleus reaction cross sections in the low energy 
range (1-100 MeV). 



1. Introduction 

Neutrino-nucleus reactions in the low-energy 
region are important for multiple reasons (see 
e.g., [1,2]). First, neutrino-nucleus reactions play 
important roles in many astrophysical processes 
including stellar nucleo-synthesis. Second, the 
observation of neutrinos emitted in various as- 
trophysical processes provides valuable informa- 
tion in both astrophysical and particle-physics 
contexts, and most of terrestrial experiments to 
measure the astrophysical neutrinos use nuclear 
targets. Third, the investigation of the neutrino 
properties with the use of accelerator neutrinos 
also very often employs nuclear targets. All these 
studies call for sufficiently accurate estimates of 
the cross sections for various neutrino-nucleus re- 
actions. It is obviously beyond the scope of this 
short note to discuss all of these individual cases. 
We instead concentrate on several selected exam- 
ples which help convey the basic features of cal- 
culations involved in low-energy neutrino-nucleus 
reactions. Since the description of a j/-nucleus 



reaction becomes increasingly complicated as the 
target mass number A gets larger, it is illuminat- 
ing to first discuss the v-d reactions {A = 2) as 
the theoretically most tractable case. We then 
proceed to give a brief discussion of the v- 12 C 
and i/- 16 reactions and finally we add a few 
comments on the neutrino reactions on medium- 
heavy and heavy nuclei. We start with the expla- 
nation of two theoretical frameworks (SNPA and 
EFT) that are of general relevance in describing 
electroweak processes in light nuclei. 



2. Standard nuclear physics approach 
(SNPA) 

In nuclear physics, the phenomenological po- 
tential picture has been highly successful. In this 
picture an A-nucleon system is described by a 
Hamiltonian of the form 



(i) 



i<3 
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where U is the kinetic energy of the i-th nu- 

rphen 



cleon, V{„ 



is a phenomenological two-body po- 
tential between the i-th and j-th nucleons, and 
the dots represent potentials involving three or 
more nucleons (which are much less important 
than V/' i h ' en ). The nuclear wave function |\Jj> is 
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obtained by solving the Shrodingcr equation 

H\V>= E\^> . (2) 

Traditionally, finding useful truncation schemes 
for | 'J' > has been an important branch of nuclear 
physics (the shell model, cluster model, etc.) but, 
thanks to the progress in numerical computation 
techniques [3] , it is now possible to directly solve 
eq.(2) for light nuclei (A<11). The phenomeno- 
logical nature of V? hen is reflected in the fact that 
its short-distance behavior is model-dependent. 
(The large-distance behavior of V? hen is con- 
trolled by the requirement that Vij should ap- 
proach the Yukawa potential.) One therefore as- 
sumes a certain functional form for V^ hen and 
adjust the parameters appearing therein so that 
the solutions of eq.(2) for the A=2 case repro- 
duce the two-nucleon observables. There are by 
now a number of so-called modern high-precision 
phenomenological N-N potential that can repro- 
duce all the existing two-nucleon data with nor- 
malized x 2 values close to 1 [4]. To describe nu- 
clear responses to external electroweak probes, we 
use transition operators that consist of the dom- 
inant one-body terms (or the impulse approxi- 
mation (IA) terms) and exchange-current (EXC) 
terms, which represent the contributions of nu- 
clear responses involving two or more nucleons. 
These transition operators are derived in such a 
manner that they are consistent with the nuclear 
Hamiltonian in eq.(l) and satisfy the low-energy 
theorems and current algebra [5]. The theoreti- 
cal framework summarized here is a cornerstone 
of contemporary nuclear physics (for a review, see 
e.g., [3]), and it is becoming common to refer to it 
as the standard nuclear physics approach (SNPA), 
the term apparently first used in [6]. SNPA has 
been scoring great phenomenological successes in 
correlating and explaining a vast variety of nu- 
clear phenomena [3]. 

3. Effective field theory (EFT) 

A new approach based on effective field the- 
ory (EFT) is rapidly gaining ground in nuclear 
physics. The basic idea is that, in describing 
phenomena characterized by a typical energy- 
momentum scale Q, we need not include in our 



Lagrangian those degrees of freedom that pertain 
to energy-momentum scales much higher than Q; 
they can be "integrated out" with only active low- 
energy degrees of freedom (or effective fields) re- 
tained. The effective Lagrangian, C e g, govern- 
ing low-energy dynamics turns out to be given 
by the sum of all possible monomials of the effec- 
tive fields and their derivatives that are consistent 
with the symmetry requirements of the original 
Lagrangian. Since a term involving n derivatives 
scales like (Q/A) n (A is a certain cutoff scale), 
the terms in £ e ff can be organized into a pertur- 
bative series in which Q/A serves as an expan- 
sion parameter. The coefficients of terms in this 
expansion scheme are called the low-energy con- 
stants (LECs). If all the LEC's up to a specified 
order n can be fixed either from theory or from 
fitting to the experimental values of the relevant 
observables, £ e ff serves as a complete (and hence 
model-independent) Lagrangian to the given or- 
der of expansion. These considerations, applied 
to hadronic systems, lead to an EFT of QCD 
known as chiral perturbation theory (%PT). A 
variant of %PT, called heavy-baryon chiral per- 
turbation theory (HBxPT) [7], is used for a sys- 
tem involving a nuclcon. However, HBxPT can- 
not be applied in a straightforward manner to 
nuclei because the existence of very low-lying 
excited states in nuclei invalidates perturbative 
treatments. Following Weinberg [8], we avoid this 
difficulty as follows. We classify Feynman dia- 
grams into two groups. Diagrams in which every 
intermediate state has at least one meson in flight 
are categorized as irreducible, and all other dia- 
grams are called reducible. We apply the chiral 
counting rules only to irreducible diagrams. The 
contribution of all the two-body irreducible dia- 
grams (up to a specified chiral order) is treated 
as an effective potential (to be denoted by V^ FT ) 
acting on nuclear wave functions. Meanwhile, the 
contributions of reducible diagrams can be incor- 
porated by solving the Schrodinger equation 

# EFT |* EFT >= £|* EFT > , (3) 
where 

A A 

H EFT = £^+£^ PT , ( 4 ) 

i i<j 
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We refer to this two-step procedure as nuclear 
xPT, or, to be more specific, nuclear xPT in the 
Weinberg scheme. 4 

To apply nuclear xPT to a process that involves 
(an) external current (s), we derive a nuclear tran- 
sition operator 7~ EFT by evaluating the complete 
set of all the irreducible diagrams (up to a given 
chiral order v) involving the relevant external cur- 
rent (s). To preserve consistency in chiral count- 
ing, the nuclear matrix element of T EFT must be 
calculated with the use of nuclear wave functions 
which are governed by nuclear interactions that 
represent all the irreducible A-nucleon diagrams 
up to v-th order. Thus, a transition matrix in 
nuclear EFT is given by 

Mf t T = < *; FT |t eft i* fft > , (5) 

where the superscript "EFT" means that the rel- 
evant quantities are obtained according to EFT 
as described above. If this program is carried out 
exactly, it would constitute an ab initio calcula- 
tion. In actual calculations, however, it is often 
very useful to adopt in eq.(5) wave functions ob- 
tained from SNPA. This eclectic approach, called 
hybrid EFT, can be used for complex nuclei (A = 

3, 4, ...) with essentially the same accuracy and 
ease as for the A=2 system [11,12]. 

4. v-d reactions 

As is well known, the v-d reactions are of par- 
ticular importance in connection with the SNO 
experiments [13]. On the theoretical side, because 
of their relative simplicity the v-d reactions serve 
as pilot cases for demonstrating the reliability of 
the available calculational frameworks. 

Nakamura et al. [14,15] performed detailed 
SNPA calculations of the cross sections for the re- 
actions: v e d — > e~pp, v x d — ► v x pn, v e d — > e + nn, 
v x — » v x pn (x=e, fi or r). The strength of the 
exchange current (dominated by the A-particle 
excitation diagram) was fixed by fitting the ex- 
perimental value of the tritium /3-decay rate . 
The results of Nakamura et al. are considered 
to be reliable at the 1 % level in the solar neu- 
trino energy range (E v < 20 MeV), and at the 5% 

4 For an alternative form of nuclear EFT based the KSW 
scheme [9], see, e.g., Rcf. [10]. 



level up to the pion-production threshold energy. 
The basis for this statement is as follows. First, 
many electroweak observables in light nuclei cal- 
culated with SNPA indicate this degree of relia- 
bility of SNPA [3,16]. Second, for the solar energy 
region, the SNPA results are supported by EFT 
calculations as well. Butler, Chen and Kong [17] 
carried out an EFT calculation of the v-d cross 
sections, using the KSW scheme [9]. Butler et 
al.'s results (after one unknown LEC is adjusted) 
are consistent with those of Nakamura et al. [15]. 
Further support comes from Ando et aL's calcu- 
lation [18]. As mentioned, hybrid xPT allows us 
to treat nuclei with different mass numbers on 
the equal footing. Park et al. [19,20], applying 
hybrid %PT to the GT transitions in the A=2, 
3 and 4 systems, noticed that only one unknown 
LEC (denoted by cIr) appears, and that da can be 
determined with good precision from . Hybrid 
XPT used in this manner is called EFT* [21,22]. 5 
Ando et al. [18] used the same EFT* approach to 
carry out a parameter-free calculation of the v-d 
cross sections, and their results agree with those 
of Nakamura et al. [15] within 1 %. Thus, at least 
in the solar neutrino energy region, the SNPA cal- 
culation has solid support from a more fundamen- 
tal approach (EFT or EFT*), and the nature of 
SNPA is such that its validity is expected to ex- 
tend up to the pion production threshold (E v < 
130 MeV). 6 

5. v- 12 C reactions and v- 16 reactions 

As the mass number A of the nuclear target 
gets larger, it becomes difficult to carry out full 
SNPA calculations of v-A reactions, which forces 
us to introduce certain approximate treatments 
of the nuclear wave functions. To illustrate issues 
involved in these approximations, we discuss here 
the v- 12 C reactions v- 16 reactions. It is to be 
noted that 12 C is an ingredient of liquid scintil- 
lators and 16 O is the basic component of water 

5 It is also called "more effective" effective field theory 
(MEEFT) [12]. 

6 The tabulation of the total and differential cross sections 
for the v-d reactions calculated by Nakamura et al. [14,15] 
is available at: 

<http://boson.physics.sc.edu/~gudkov/NU-D-NSGK>, 
and <http://www-nuclth.phys.sci.osaka-u.ac.jp/top>. 
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Cerenkov counters. 

The best studied method for truncating nuclear 
wave functions is the shell model, in which we in- 
troduce an average single-particle potential V ave 
to generate single-particle orbitals from which A- 
body wave functions are formed. If we use as V ave 
the harmonic oscillator (HO) potential (known to 
be a good approximation for describing bound 
states in light nuclei) , the degree of truncation of 
nuclear wave functions can be specified by stat- 
ing up to what excited HO configurations (nhcu 
excitations, where u) is the HO angular frequency, 
and n — 0,1,2, ...) are included in the calcula- 
tion. With the use of a truncated model space, 
the nucleon-nucleon interactions as well as tran- 
sition operators for various nuclear observables 
get renormalized into effective operators, which 
may be schematically denoted by V^/J and T e f$ , 
respectively. (This is similar to what happens 
with the introduction of an EFT.) There exist 
two ways to get information on V%fjq and Tiff- 
in one method they are derived from the "bare" 
forms (pertaining to untruncated cases) by incor- 
porating the contributions of virtual states which 
lie outside the chosen model space [23]. In the 
second method, we simply list up all formally al- 
lowed terms for or %, ff and determine their 
prefactors by fitting to the available experimen- 
tal data. Insofar as the number of parameters to 
be fixed is sufficiently smaller than the number 
of data points, this empirical effective operator 
method (EEOM) has predictive power [24-26] . 

For closed-shell light nuclei such as 16 0, one 
can handle a relatively large model space; for in- 
stance, a shell-model calculation including up to 
47iw configurations was reported [27]. For nu- 
clei away from shell closures the situation is less 
favorable, and one is forced to use certain ap- 
proximations instead of full shell-model diagonal- 
ization. 12 C belongs to this category, and the 
random phase approximation (RPA) is commonly 
used for this case [28]. (See however a recent no- 
core shell model calculation for the A=12 system 
by Hayes et al. [29].) 

In assessing the reliability of using the shell 
model and its approximations in calculating v- 
nucleus reactions, we need to pay attention to an 



apparent anomaly known for the inclusive reac- 
tion 12 C(^,^) 12 N* [30]. According to an LSND 
experiment [30], for 123.7 MeV<^<280 McV, 
the measured flux-averaged inclusive cross sec- 
tion was more than a factor of 2 lower than that 
predicted by the Fermi-gas model [31] 7 and by 
an RPA calculation [35]. Although the incident 
neutrino energies in this experiment lie above the 
energy region mentioned in the title of this talk, 
the range of nuclear excitations involved in this 
experiment are relevant to lower-energy (E v < 100 
MeV) neutrino reactions as well. It therefore con- 
cerns us whether there is really an anomaly in the 
12 C(f /i , /^) 12 N* reaction. When the final report 
on the LSND experiment appeared [36] , the mea- 
sured flux-averaged cross section was still signifi- 
cantly lower than the existing theoretical predic- 
tions of the time. Most recently, however, there 
has been much progress on the theoretical side. 
Within RPA, Krmpotic et al. [37] report that 
the implementation of particle-number conserva- 
tion in quasi-particle RPA brings the theoreti- 
cal prediction into agreement with experiment. 
Meanwhile, Oset and his collaborators have de- 
veloped a formalism which incorporates long- 
range nuclear correlations via RPA, and which 
in addition accounts for final-state interactions 
and the Coulomb corrections. Their very recent 
work [34] reports that there is no anomaly in 
the C(v^, /i) 12 N* reaction. Furthermore, in a 
no-core shell model calculation [29], Hayes et al. 
point out that the inclusion of a realistic three- 
body interaction can remove the discrepancy be- 
tween theory and experiment. All these develop- 
ments indicate that the long-standing "anomaly" 
is finally gone and that there is no basic prob- 
lem with the currently available nuclear physics 
approaches. 

Haxton [38] was the first to point out the im- 
portance of the is- 16 reaction in water Cerenkov 
counters, and he calculated the relevant cross sec- 
tion, using 2huj shell-model wave functions and 
the transition operators whose strength are renor- 
malized according to EEOM. Here, the inclu- 
sive nature of the reaction under consideration is 

7 A recent study [32] has questioned the quantitative relia- 
bility of the Fermi-gas even in the quasi-elastic scattering 
regime. See also [33,34]. 
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taken into account by summing all energetically 
allowed final bound states of the A=16 nuclear 
system. Interestingly, Kuramoto et al.'s calcula- 
tion [39] based on the relativistic Fermi gas model 
(RFGM) [40] indicates that the RFGM cross sec- 
tions and Haxton's results can be smoothly ex- 
trapolated into each other at around E v = 60 
MeV; see also Bugaev et al. [41]. The most recent 
examples of a calculation for the i/- le O reaction 
can be found in [42,34]. 

6. Neutrino reactions on medium-heavy 
and heavy nuclei 

The study of low-energy neutrino reactions on 
medium-heavy and heavy nuclei is of great cur- 
rent importance in various contexts. For one 
thing, some nuclei of this category are either used 
or planned to be used in solar neutrino exper- 
iments with various specific features. Further- 
more, many reactions of this type play important 
roles in supernova explosions and neutrino nucle- 
osynthesis, see e.g. [43]. It is to be noted that, de- 
pending on the condition of a core that occurs in 
stellar collapse, the average electron-neutrino en- 
ergy can range typically from 10 to 50 MeV, and 
that cross sections for v e capture on iron-group 
nuclei through A=100 are needed to accurately 
simulate core depleptonization and to accurately 
determine the post-bounce initial condition [43]. 

Providing sufficiently accurate estimates of the 
cross sections for all the relevant neutrino-nucleus 
reactions is one of the imminent challenges facing 
nuclear physics. To push ultra-large scale shell- 
model calculations (or, if possible, even SNPA 
calculations) for higher and higher mass numbers 
is an important direction of theoretical efforts. 
However there will be quite sometime before these 
frontal attacks can cover all the relevant nuclides. 
The best strategy one could take at present is to 
use a model (such as RPA or EEOM) that can 
cover a wide range of nuclides and to improve or 
gauge its reliability with the help of experimental 
data. For instance, Langanke et al. [44] consid- 
ered the possibility of utilizing EEOM and high- 
resolution electron scattering experiments to cal- 
ibrate neutrino-nucleus cross sections relevant to 
supernova neutrinos. It is also noteworthy that 



there exists an experimental project [45] which, 
taking advantage of ultra-high-intensity neutrino 
beams that will become available at the Spalla- 
tion Neutron Source (SNS) at the Oak Ridge Na- 
tional Laboratory, aims at measuring the cross 
sections of several neutrino- nucleus reactions with 
hitherto unachievable precision. These experi- 
ments are expected to shed much light on the reli- 
ability and the predictive power of nuclear models 
to be used for studying supernova explosion and 
nucleosynthesis. 
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